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Abstract 

Starting from a microscopic approach, we develop a covariant formalism to describe a set 
of interacting gases. For that purpose, we model the collision term entering the Boltzmann 
equation for a class of interactions and then integrate this equation to obtain an effective 
■ macroscopic description. This formalism will be useful to study the cosmic microwave back- 

ground non-perturbatively in inhomogeneous cosmologies. It should also be useful for the 
study of the dynamics of the early universe and can be applied, if one considers fluids of 
galaxies, to the study of structure formation. 
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1 Introduction 

i 

^Jy The formation and evolution of cosmological perturbations, leading for instance to the cosmic 

microwave background fluctuations, are usually studied in the framework of cosmological pertur- 
bations about a Friedmann-Lemaitre spacetime In such a framework, the matter is described 
by a set of decoupled species including cold dark matter, neutrinos and a baryon-photon fluid 
(since one must consider the coupling between these two species via Compton scattering). In 
the latter case, the interaction between the two coupled components (i.e. photons and baryons) 
is studied via the linearised Boltzmann equation. A first generalisation that comes to mind is 
to develop a framework to model this interaction in a covariant way, i.e. without specifying 
the spacetime geometry. Such an approach will be useful for a study of the cosmic microwave 
background in inhomogeneous cosmologies [Q, [§ • Earlier in the history of the universe all these 
fluids were coupled through scattering and/or annihilation process Q. This leads to the second 
generalisation which is to model this coupling in a very general way for a wide class of interac- 
tions. Such a formalism will then be adapted to study the dynamics during the early universe, 
independently of the spacetime geometry. Moreover, it can be applied to structure formation, 
since it can describe fluids of galaxies. 

The Einstein equations and the Bianchi identities lead to the fact that the total energy- 
momentum tensor of matter is conserved 

V M T^ = 0. (1) 

When the matter content of the universe is a single perfect fluid, this equation splits into two 
equations, the Euler equation and the matter conservation equation. The latter can be reset as 

V M n^ = 0, (2) 
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meaning that the number flux vector, n^, is conserved. 

If one now considers a system composed of an arbitrary number of fluids, this property applies 
to the "global fluid" and to each fluid only if they are interacting by gravitational interaction 
(i.e. if they are decoupled). If we want to study the evolution of one particular fluid of a system 
of interacting fluids, we have to take into account a force and the preceeding law can be rewritten 
as 

V M lf = F?. (3) 

Furthermore, if the collisions are not elastic, then the total number of each species will not be 
conserved and we will have to take into account a source term 

VX = <*■ ( 4 ) 

Such a splitting of the global energy-momentum tensor (T^ u = ) and of the number 

flux vector (n^ = J2nf) is not straightforward. In fact, for a general multi-constituent fluid || 
the lagrangian of the system, C say, will not split in the sum of individual lagrangians describing 
the dynamics of each fluid. Since the energy-momentum tensor is related to this lagrangian via 

TIXV= 2 8£J=a 



-g $g, 



one cannot define the energy-momentum tensor of a single component of the fluid system. 
However, as long as we are dealing with gases, such a splitting will be possible. Indeed, each 
fluid is then conveniently described by its distribution function, /j say, from which one can define 
an energy-momentum tensor (see &|| and &j3| for the complete definitions) as 

if = /pWi(x,p)7r+(p). 

This standard kinetic theory approach H will be valid if the particles interact weakly, i.e. for 
systems which are not too dense. As explained in the section 5.4, such a description is adapted 
to cosmology. 

A way to model the interacting forces, Fj 1 , and the source terms, e^, is to relate them to the 
collision term entering the equation of evolution of the distribution function. We will then start 
from the kinetic theory, try to give a general form of the collision term for a class of interactions, 
compute the forces and the source terms, and end up with a macroscopic description of the 
dynamics of interacting fluids. It will be non perturbative, in the sense that we will not linearise 
the metric, and will not assume anything on the form of the fluid (e.g. perfect fluid). 

Both the kinetic theory and the theory of fluid dynamics have been studied in the context of 
general relativity. On the one hand, a general covariant formulation of kinetic theory in general 
relativity was first developped by Tauber and Weinberg |7| and independentely by Chernikov 
||] (see also Marie || and Israel and Stewart flOfl ). These authors were mostly concerned 
with relativistic generalisations of classical gas theory (proof of the H-theorem, equilibrium 
configurations etc.). Many authors have also studied the coupled Einstein-Boltzmann equation 
and solved it |ll], |l2| in some cases. It is not our goal to solve such a general problem here, 
since we do not want to assume anything concerning the symmetries and/or the geometry of 
the spacetime. 

On the other hand, relativistic fluid dynamics has been studied in detail by many authors 
(see e.g. |l3|), who usually do not deal with a system of interacting fluids even if they 
present a general formalism for a multi-constituent fluid. In fact, as will be explained later, the 
formalism developped by Carter || for multi-constituent perfect fluids does not overlap with 
the description which we shall give here. 
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In this article, we stand in between the kinetic and the fluid descriptions. Note that even 
if such attempts have already been made, they relied on a different approach to the problem. 
For instance, Lindquist 14] studied the diffusion of photons under the transport approximation 



(i.e. he studied the evolution of particles flowing through an emitting and absorbing medium 
described in macroscopic terms). Most of the studies of the transfer equation are also based 



on moment methods ||, on the "grad's method of moment" [15 or on a spherical harmonic 



analysis [16]. We will not use such expansion methods here. 

We will first use the microscopic approach (& ||) to relate the interacting force and the source 
term to the collision term of the Boltzmann equation (& |3|). In & [|, this result is used to compute 
the force on a system of conducting fluids. We then turn to the computation of the force for 
elastic collisions in & |j[ where we discuss the general form of the force and then compute it for 
cases of cosmological interst (e.g. photon-fermion (using results on the photon-electron collision 
|jj"7|| ) and fermion-fermion scattering) . & ^ is devoted to the computation of the general form of 
the source term and of the force for inelastic collisions, namely fusion and fission, and we use it 
in the case of photon (bremstrahlung) of fermion-antifermion annihilation and of recombination 
(which is very important for a study of decoupling in cosmology). 

This general formalism is then used in to establish the equation of evolution of a Compton 
scattering coupled photon- matter system in Cosmology to linear order (& |7|) . All the theory of 



multifluid linear cosmological perturbations [18] can be recovered from our formalism which is 
non perturbative. 

2 Definition - notations - microscopic quantities 

The goal of this section is to introduce the distribution function and its equation of evolution. 
For that purpose we need to introduce the space on which such a function is defined. We finish 
by the description of the 3+1 splitting with respect to an arbitrary vector field. 

2.1 Distribution function 

Let us first consider a single test particle with mass m which moves in a gravitational field. Its 
motion is determined by the geodesic equation 



dx^_ 1^ ^ 
P ~ dX' dX dX +L -P pp " U ' w 



where A is an affine parameter defined by the requirement that p^ 1 be the 4-momentum. Here- 
after, V„ denotes the covariant derivative associated to the metric <7„„, whose Christoffel symbols 
are T% p . Let us note that if there are non gravitational forces (e.g. electromagnetic forces) then 
we have to modify this equation (see & |j) . 
The rest mass of the particle is defined as 

m? = —p^Pu- (6) 

Thus, according to (|5|), the state of the particle is determined by the couple (a^p^) and the 
phase space is then the tangent bundle over the spacetime manifold, i.e. 

T = {(x,p),xGM, P GTJ, (7) 

where A4 is the space-time and T x is the tangent space to M. at x. From now on, we will use the 
bold style to denote quadri-vectors when components notations are not needed (hence p = p^ L 
and p 2 = p^p^ and greek indices run from to 3). 
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The volume element on T x supported by the displacements dpi,dp2,dp3,dp4 (with compo- 
nents dp" etc.) is 

tt(p) = e a {3 lS dpfdp%dpldp 5 4 , (8) 

where e\ a p^ is the totally antisymmetric tensor such that £0123 = \J— 9- 

We also define 7r + (p), the volume element corresponding to the subspace of T x such that p^ 
is non-spacelike and future directed, 

7T + (p) = H(-p^)H (-p 2 Mp), (9) 

where H is the heavyside function (i.e. H(x) = if x < and H(x) = 1 if x > 0) and an 
arbitrary timelike vector field. 

T x is sliced in hypersurfaces, ? m , of constant m called the mass-shell, and defined by 

P m (x) = {p G 7i,p% = -m 2 ,p% > o} . (10) 

The volume element (g) on T can then be decomposed on a volume element, rmr m , on P m by 

7r + (p) = mir m (p)dm. (11) 

The factor m allows one to include particles of zero rest mass (see Ehlers ||). This defines the 
induced volume element m7r m (p) on V m . 

If we introduce an arbitrary future directed unit timelike vector (i.e. satifying u^u 11 = — 1), 
the 3- volume supported by the three displacements dx%, dx2,dx^ (with components cfccf etc.) in 
the hypersurface perpendicular to is 

dV(u) = exa^^dxidx^dxl- (12) 

We now consider a single fluid composed of particles which mass are a priori different. The 
distribution function, /(x, p) will be defined as the mean number of particles (on a statistical 
set) in a volume dV around x and 7r(p) around p measured by an observer with 4-velocity u^, 

dN(x, P ) = /(x, p)(-jAi M )dy(u)7r(p). (13) 

The assumptions involved in its existence have been discussed in details by Ehlers . Synge p0| ] 
has demonstrated that (— p^u^)dV{u) is independent of u^, which implies that the distribution 
function is a scalar. Moreover, /(x, p) > for all x^ and all allowed p^. 

For a gas, dN is the number of particles in a volume dVir(p) thus the smoothness of / 
depends on the existence of a sufficient number of particles. 

2.2 Equation of evolution 

The equations of motion (|5|) define on T an operator called the Liouville operator (see e.g. |2^] ) 
which reads 

p dx^ dX dp* dA' 1 ' 

which characterises the rate of change of / along the particle worldlines. Using (||), this operator 
can be rewritten as 

C[f]=p^d,f-T^p^f. (15) 

The fact that the mass m of the particle as defined in (|6|) is a scalar function which is constant 
on each phase orbit leads to 

C[m 2 } = 0. (16) 
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The Boltzmann equation states that this rate of change is equal to the rate of change due to 
collisions, i.e. that 

£[f] = C[f]. (17) 

C[f] is the collision term and encodes the information about the interactions between the par- 
ticles of the fluid. 



If we know consider a system of N fluids (labelled by each of which is described by 

its distribution function /j(x, p), the Boltzmann equation for a given fluid i becomes 

nn\ •/•-./): cii- (is) 

3 

Cj[fi, fj] is the collision term describing the interaction between the fluid i and the fluid j. For 
elastic collisions, it must satisfy the symmetry 

r ; /„./• ; r, /,,./• ;. (19) 

which means that in a collision between i and j the two distribution functions undergo the same 
change. Following Israel and Stewart |l(J, we will require that Cj[fi,fj] is a local function of 
the "/j" (i.e independent of their derivatives). 



2.3 3+1 splitting 

We perform a 3+1 splitting with respect to an arbitrary timelike unit vector field (we discuss 
which vectors to use according to the problem at hand in & |3.2[ ). The projection tensor into the 
"rest-space" of an observer moving with this 4-velocity is defined by 

-L/^= g^u + U^Uy. (20) 

Any vector p^ can be decomposed with respect to as 

= Xe^ + Eu" , (21) 

where the energy (E), the norm of the particles' 3-momentum (A) and the direction with respect 
to are 

£7 = (-p%), A 2 El F // = J? 2 -m 2 and ^ v P " = , (22) 

so that e M and satisfy 

u% = 0, = 1. (23) 

We have dropped the u-dependence since there is no ambiguity. However we will restore it if 
necessary (e.g we will write -L^f). Here, we have split the 4-momentum but any other vector 
can be split according to the same procedure (e.g. in equation 38). Note that for a zero- mass 
particle these relations lead to A = E. 

From these definitions and equations (|9|) and (0), we can reexpress the volume elements on 
T x and on the shell-mass as (see ffl, [l6| for details) 

vr + (p) = mdmXdEdil <==^- 7r m (p) = XdEdft, (24) 

where dQ is the solid angle spanned by three independent e^. The passage from (H) to (p3) can 
be seen as the passage from a cartesian coordinate system to a spherical coordinate system (i.e. 
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we replaced the integrations on dp a by an integration on the norm and the angles.) 
We also introduce the symmetric traceless tensor 



which verifies 



A^J_^=0, A M ^ = 0, A liv u" = and A^e" = -e v . (26) 



(27) 



To finish, we give the following useful integrals (see e.g. |l( 

f w = / n = 2p + l 

J 6 '" 6 " 4tt \ jjjj (l^w ... n = 2p, 



where satisfies equation (p3|). In particular we have 



4, ^ X ' < 28 > 

f, = 5^. < 29 > 

f A a/3 A 7^ _ J_ C_ 2 _L"/3_L75 +3 _|_a7_|_/3<5 +3 _|_orf_|_a</\ _ (30) 

J 4ir 45 V / 

We also choose the following conventions of symmetrisation and anti-symmetrisation 

A (ai...a„) = — { A a(i)...a(n) and A [ai...a n ] = H e ( a ) A a(l)...a(n) > 

crGperm(l..n) <rGperm(l..n) 

where e(er) is the signature of the permutation. 

3 Macroscopic quantities 

In the previous paragraph, we have given a microscopic description of a set of interacting gases. 
The goal of this section is to define a set of macroscopic quantities from the distribution function 
and the collision term and then find the relations between these quantities. 

3.1 Definition 

At any point x, one can introduce, following Ellis et al. given a distribution function fa a 
set of macroscopic quantities associated with each fluid i by 

X^-^(x) = / (-p^)°/V 1 -P Mn /i(x,p)7r + (p) = / / mV 1 ...^/ I (x,p)m^ ( i m 

JT X Jm J Vm 

(31) 

and 

^f„ 1 - Mn (x)= / (-^) a/2 ^ 1 ...p w "C[/i](x,p)7r + (p)= / / my i ../»C[/ 1 ](x,p)mf m d m , 

(32) 

where m is the mass of the particles (defined in (^)) and a an integer. The particles of a given 
fluid can have different rest mass (this is the case e.g. when one is dealing with a fluid of stars 
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or of galaxies) . We assume that each distribution function vanishes at infinity on the mass shell 
quickly enough so that all these integrals converge. 

These quantities are clearly totally symmetric and follow (from equation @) the general 
recursion relation 

g^x^-^W = -x«- 2 (x). (33) 

The functions Y follow the same recursion relation. 

The distribution function and the collision term are related via the equation ( |i~7l) which 
implies the two sets of quantities X and Y are not independent. It can be shown (see Ehlers 
H) by using the Boltzmann equation fll~7|) , (|l6| ) and that d^ir = (In ^—g) 7r(p) = T v ■n(jpi) 
that the quantities X and Y are related by 

V l *X?£"*" = Y£" ftn , (34) 
on which our following construction is based. 

Among all these quantities, some are important in many applications, 

n» = Xg, N» = XI and = X$ v . (35) 

The vector is the number flux vector which is used to define the average number flux velocity 
vector and the proper density n measured by an observer comoving with the fluid by 

n^ = nv^, v fM v fl = -l. (36) 

Similarly, is the average mass flux vector which defines the mass flux velocity vector and 
the energy density p by || 

N» = pW, V% = -1. (37) 
For an observer with an arbitrary 4- velocity u^, can be split as 

n^ = nu^+f; n = _( n ^M). f=±*f. (38) 

is the diffusion current as measured by this observer and vanishes when = v^. 

Tu V is the energy-momentum tensor. It also determines a unique average velocity, namely, 
its timelike eigenvector (see e.g. ||]). Whatever the timelike unit vector field u^, chosen as time 
direction, we can split the energy-momentum tensor under the general form 

T/u, = pu^Uy + P +2q {p u v) + vr^, (39) 

the quantities p, P, q^ and tt^ being defined as 

p = T^u u ; P = q^ = - ±^T ua u a ; ir^ =±«±? T a p - P ±„ u . (40) 

This decomposition is the most general splitting with respect to the arbitrary vector field vJ 1 of 
a tensor of rank 2. The two quantities q^ and tt^ u are respectively called the energy flux and 
the anisotropic stress and verify (from 40) 

= K = n ^ utl = °- ( 41 ) 

By using the definition of T^ v from the distribution function (equation |3l| ) and performing the 
splitting (|2l|), a simple identification with equation (|39|) easily shows that 

poo poo P 

p(x) = / dm dE f(x,p)E 2 \dn, (42) 
Jo Jm Jn 

1 poo poo p 

P(x) = - / dm dE /(x,p)A 3 dn, (43) 
3 jo Jm Jn 

poo poo p 

^(x) = / dm dE /(x,p)PA 2 e^ (44) 
Jo Jm Jn 

poo poo p 

tt^(x) = / dm dE f(x,p)\ 3 A fl „dn. (45) 
Jo Jm Jn 
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All these quantities depend intrinsically on the choice of the vector field u^ 1 via the splitting 



defined in equation (21). Thus, p, P, and ir^ will respectively be the energy density, the 
pressure, the energy flux and the anisotropic stress measured by an observer comoving with v?. 

3.2 Macroscopic fluid dynamics 

Using the equation (|34|) as well as the definitions (j3l[-|32"l) and (|35|) , we can relate the force Fj 1 
and the source term e{ to the collision term defined in (||-(|) by 



and 



if(x) = ^ =/ ^Q[/ i ](x, P h + (p), (46) 

e, (x) =Y it0 = f d [M (x, p)vr + (p) . (47) 

The Bianchi identies state that the total energy-monentum tensor is conserved, which implies 
that we must have the usual action-reaction law, i.e. that 

E F ? = 0, (48) 

i 

as long as there is no long range external force (such a force can only be of electromagnetic 
origin |J ; this will be studied in & ||) . 

Most of the time we will have to pick up a special frame to compute the collision term. Some 
choices are possible even if there are not compulsory. 

• When a massless particle is interacting with a massive particle we will choose the rest 
frame of the massive particle and thus 

vP = —. (49) 
m 

• When one has an elastic collision of two massive particles we can use the center of mass 
rest frame defined by 

P» =rf+ pit and = ■ (50) 

Besides these two velocities, we have seen that there exist some preferred timelike vector fields 
associated with the motion of the matter (e.g. v^ 1 , V^). 

One could add to this dynamical description a thermodynamical description. We have to 
emphasize here that we must give an equation of state to close the system. This comes from 
the fact that when developped into moments, at a given order the Boltzmann equation involves 
multipoles of higher orders |15]| . To close the system one has either to troncate the system at 
a given order (with all the arbitrariness it implies) or give an (or more if needed) equation of 
state for the fluid. Even if we do not specify it, we assume that such an equation can be given 
for concrete applications (see & 0) . 

The goal of the following section is to compute explicitely the quantities (E^-47) in terms of 



the macroscopic variables defined in section 3.1 
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4 Conducting fluid 



The easiest case where one can compute the force acting on a fluid is the case of a conducting 
fluid in a electromagnetic field. This has been computed for a single fluid by many authors (see 
[13] for a review concerning electrodynamics in continuum media). We will just make the link 
between the microscopic and macroscopic approaches and show how useful the latter can be. 
Let us start from the usual approach with a single fluid. 

Since the total energy-momentum tensor, TJ^, is conserved, we have 

VJTV = -V U TZ, (51) 
the electromagnetic field energy-momentum tensor being defined by 

TZ = T^n-\g^T\ (52) 

where is the electromagnetic tensor. Using Maxwell's equations, 

= j», (53) 

j* 1 being the current density, the electromagnetic force on the fluid is 

Fj> = -V v T% = r£j a . (54) 

However in the case of a multi-fluid system, the only force that can be computed with this 
method is the global electromagnetic force on the system of fluids and not the force on each fluid. 

If we turn to the microscopic approach, then we have to take into account the fact that, 
because of the Lorentz force, the particles do not follow a geodesic between two collisions. Their 
equation of motion is then given by 



tfVX = e^p u <j=> + T^pP = e^p u , (55) 



dp' 
~Jt 

where ej is the charge of the particles. 

If we cast this relation in the Boltzmann equation as we did in the section 2.1, it can be 
written 

£[/] = C B [f] = - ei p^^. (56) 
Hence the electromagnetic force acts as a collision term. In the previous equation C stands for 



the Liouville operator with no electromagnetic field, as defined in section 2.1 
We can now compute the force coming from this collision term, 



df_ 

dp?' 

e^nf^FZj?, (57) 



where we have performed an integration by part. Thus, the force acting on a single fluid can be 
expressed in term of macroscopic quantities, namely the electromagnetic tensor and the current 
density. 

If we now compute the total electromagnetic force on all the fluids we get 

F B=J2 F B-*i = ^a^Ji = F%J a , 
i i 

which is the result we obtain with the macroscopic approach. 
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5 Elastic collisions 



In this paragraph, we will compute the force for elastic collisions. We begin by a derivation 
for binary collisions of classical (in the sense that they are non quantum) particles, trying to stay 
as general as possible. We then turn to the case of the Compton scattering. In the Thomson 
limit, even if electrons and protons are not classical particles, the computation of the force turns 
out to be an application of the general case. However, in this case, we can give the form of the 
corrections coming from the quantum statistics. 

We finish by a discussion on the range of validity of these computations and determine the 
domain of applicability (range of temperature) in which they can be used and then discuss the 
pertinence of this formalism for cosmology. 



5.1 general case 

We will give the most general form for the collision term for the process 

A(p)+B(d)^A(p')+B(^). 

Before computing the force, we will study its general symmetries, which starts by the study of 
the kinematics of such a collision. During the computation, we will have to break this symmetry 
by working in the rest frame of one of the particles, and restore it at the end. 

We first work in the rest frame of the center of mass and thus use the vector field defined 
in ( |50|) to perform the splitting. We then have from equation ( p2|) 

E A = TO > " and E B = m * ~ p P ^ with P = (m\ + m% - 2p^)V\ (58) 

A 2 ^Ai = A|= ( ^ )2 ~ 2 ^ m| , and (-p^) = E A E B - \ A X B . (59) 

E A and E B are the energies of the two particles A and B in the center of mass rest frame, A is 
the amplitude of the particles' 3-momentum is this frame. 

The vectors and e B , the direction of the ingoing particles, are given by 



= ^{E B jf-E Aq ^}. (60) 



The general form for the collision term for binary collisions of uncharged classical particle is 

§ 

C[f A )(p) = J {-fq v ) [MpOMqO - /U(p)Mq)] W(p,q,p',q / )^ + (q)^ + (q> + (p / ), (61) 

where VF(p, q, p', q') is the probability of a collision (p, q) — * (p',q'). General expressions for 
processes including electromagnetic effects and quantum statistical effects have been proposed 



||, [1I|. We will take such effects into account in the next section (see equation (|9q)). Here, we 
want to be as general as possible and try not to describe the interaction in detail but give the 
general form of the force. The microscopic reversibility of the collision imposes the symmetry 

W(p,q,p',q') =^(p / ,q / ,p,q). (62) 
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Before we estimate W(p, q, p', q'), let us study the symmetry of the force 

F^A = J (-P V Qv) {/A(pO/B(qO " /A(p)/ B (q)}p^(p,q,p , ,q / )^+(q)^+(q , )vr + (p / )^+(p)- 

(63) 

If we sum the two forces (i.e. F B ^ A and F^-*b) an d se ^ ^(p> 1; P j = 5^(P— P')i?(p, q, p', q'), 
where R and have the same symmetry, in order to make the conservation of energy-momentum 
explicit ,then since P 2 = P' 2 implies that (-pq) = (-p'q'), the integrand of F^ A + F^ B 
will be antisymmetric in the transformation (p, q) — ► (p', q') and thus 

F^ A + F^ B = 0. (64) 

This was expected by construction from ( |48|) but had to be checked on the general form. 



To compute the force, we must relate the W(p, q, p', q') probability to the differential cross 
section d& pq ^ p q defined as 

^(p,q,p',q')vr + (q')vr + (p / ) = (-^da™^' . (65) 

The cross section can be decomposed p3[ as the product of a matrix element M kp ^ k ' p ' and of 
a two body phase space element as 

1 1 



da pq ^ p ' q ' = -- \ M k P ->k'p'\2 f D k P ^k'p') ( 66 ) 



with 



(D kp - k ' p 'f = (2^)V 4 )(p + q - p' - q> + (q>+(p'). (67) 

If we consider massive particles with respective rest mass uia and mg, one can convince 
oneself that the matrix element has to be of the form 

\M kp ^ k ' p '\ 2 = aT(m 2 A ,m 2 B ,p^,e^e'n, (68) 

where a is the scalar cross section and the function T can be decomposed as 

T = a + /3e^ + 7 A^A^ + ..., (69) 

where the coeficients a, j3, 7 depend on m A , m 2 B and p v q^ and where and e'^ are the directions 
of the ingoing and outgoing particles. This form comes from a multipole expansion of the matrix 
element in which the coefficients depend only on the scalar invariants of the collision (an example 
of such a function T is given in equation (|B6|)). As will be shown below, these are the only relevant 
terms to compute the force since higher multipoles will not contribute. 

We will now assume that the particles are either non-relativistic in the center of mass rest 
frame or that one of the two particles is massless. The first approximation can be stated by 
assuming that the 3-momentum (A defined in (|59|)) is small compared to the total energy (Ea 
and/or Eb) which is of same order of magnitude that the rest mass (mA and/or mg), i.e. that 

A <C E ~ m. 

Since T(m 2 A ,m 2 B ,p IJ, q^) = T(m\, m 2 B , A, Ea, Eb), in this approximation we can Taylor expand 
each coefficient of equation (p9|) in power of X/ttiab, with vtiab = inf(mA,TnB) as (on the 
example of a) 

a(m 2 A ,m 2 B ,p^q^) = a NR ^ — a + O ' 



m.AB \ \ m AB 
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On the other hand if we assume that one of the particle is ultra-relativistic, then we can make 
a "Thomson-like" approximation, which says that the energy of the zero-mass particle is small 
compared with the rest mass of the particle it scatters with (i.e. p^q^ <C rn 2 AB with ttlab = 
sup(rriA,mB))- Each coefficient of equation ( |69|) can then be expanded as 



a(m A ,m B ,p 1 q il ) = a UR ^f-a + CM — ^ 

m AB V V m AB , 

Technically these two approximations reduce the knowledge of the function T to a set of scalars 
( (&nr, Pnr, Inr, oe, ft, 7) in the first approximation and (auR> Pur,7UR,&> m the second 
one). 

Since at zeroth order the two approximations lead to the same form of the force, we will first 
compute this term and then evaluate the first order correction in the two approximations. Note 
that even if the form of the force is the same at zeroth order, the force itself will be different 
because a^R and auR are not the same constants. We have 

onr = a(m 2 A ,m B ,E A = m 2 A ,E B = m B , A = 0), 

auR = a(m A ,m 2 B ,p^q^ = 0). 

For the sake of simplicity we set a = (chnr, our)- 

We will compute F AB in the rest frame of B. As explained at the beginning of this section, 
this breaks the symmetry between A and B (see (p^)). Therefore, we will have to restore this 
"hidden" or "lost" symmetry at the end of the computation. Another solution would have been 
to split the force into two halves and compute one part in the rest frame of A and the other 
half in the rest frame of B. Let us note that this symmetry can be restored only in the first 
approximation where we can compute the force either in the rest frame of A or the rest frame 
of B. In the second approximation, this is no longer possible and the force will not be obviously 
symmetric. We define F AB by 



AB 



HQ' 

(~p u Qu) {/a(p') - /a(p)} /fl(qKT(mi,m|,p^ M ,^)— 7r+(q)AdEdfi 



rfB 

(70) 

where ll rfB" means that the quantities are evaluated in the rest frame of B. From equation 



(64), it follows that in the first approximation (i.e. the non-relativistic approximation) we have 



r A^B — r [AB]i 

(where [AB] means that we anti-symmetrise on A and B the expression (fTOj) ) and in the second 
approximation (the ultra-relativistic approximation) 

^a^b = Fabi ^b^a = ~^AB- 



We perform the splitting (|2l| ) of p^ with respect to defined by = and inject in the 
integral (|70|). For the sake of clarity we will split F AB in 

F AB + F AB where (°) and (1) refer 
to the zeroth and first orders in the expansion of T either in "A/mAB" or in "(— pq)/m AB v . 
Thus, taking into account the fact that q^e^ = and that / e^dQ = J A^dQ = 0, it follows 

Fab = aa J PVWA(p)/B(q)7r + (qh+(p) 
+ a J (SV + E\m B e») / B (q)7r + (q) 

{« / fAW^+Wj WKa^+7A^ / f A (p')A' ap ^XdEdn. (71) 
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Using the definitions (4^-45) and the integrals (28-3C), this reduces to 

F { Z = a \an Bv T% + ap A n B + 3^^} ■ (72) 

Let us recall that q A is the energy flux with respect to n B , i.e with respect to the unit vector 
u B colinear to n^, 

QA = ~ y-A U Bu + J- A u BaUB/3UB] ■ 

Thus, we have obtained the expression of the force on the fluid composed of particles A in 
terms of the macroscopic quantities describing the two fluids and of the coefficients (a, (3, 7, a) 
describing the collision. 

Let us now turn to the evaluation of the first order corrections. We begin by the "non- 
relativistic" approximation, 

F^B = —5 /^pV/A(p)/B(q)A7r + (q)7r+(p) + — / (#V + £Am fl e^) / fl (q)7r + (q) 
mAB J rn AB J v > 

\a J fA (p>)*L + far J fA ( P y a ^ + ^J f A (p')A> a ^} x'dEdn. (73) 

Using the definition of A ( p^) , and the set of new quantities defined by 

J A = / / f A (x,p)\ 2 E 2 dEdn, (74) 
Jm Jn 

poo p 

J% = f A (x,p)\ 3 Ee a dEdn, (75) 

Jm J f2 

Jf = |p>Ve 7 /A(x,p)^ + (p), (76) 



and the definitions (42-45) and integrals (28-29), we obtain 



Fab = — {uub^a + &J A n B + fejf) • (77) 
mAB y - 



We will now turn to the second approximation (i.e. p^q^ <C m AB) anc ^ evaluate the first 
order contribution. It reads 



F ( 2b = I P fl P U P X QuqxfA(p)f B (ci)7T + (q)7T + (p) f (m B E\» + E?\m%eP) /s(q)7r + (q) 

m AB j m AB J 

{ & J f A (p')^ + Pe a J f A {p')e'J^ + 7A^/ f A (p')K^} ™EdQ, (78) 



and thus 



Fab = — — &N^I A + f - ± B a T^I Aa - ~aT BvX X^ , (79) 
m A mB 3 h 



mAm B 

where we have introduced the two macroscopic quantities 



poo p poo p 

I A = I f A \E 3 dEdtt, I A = / f A X 2 E 2 e a dEdQ, (80) 

Jm Jfl Jm JO. 



and where the indice B in A. B a means that we used the 4- velocity u B . X^ 1 ^ is defined in (31). 
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In conclusion the force is given by 



B^A ~ I \ t [AB] +t [AB}) 

F^A = (F^ + F^), (81) 

according to the first or the second approximation. is given by © an F* 1 ) either by @ 
or ©. 

Let us stress that we could go on in the expansion of T and compute corrections to the force 
at different orders. This would however involve the introduction of new macroscopic quantities 
(like e.g. J a, J%---) an d of tensors of higher rank (like e.g. Xq^ X ...). This is an example of what 
we mentioned in & 13.21 since F<® involves <f and F^\ X^ v . 

In the case of quantum particles, we have to take into account quantum statistics effects 
which will be evaluated in the next section. 



5.2 Photon-electron scattering 

We will focus here on the elastic Compton scattering between electrons and photons, 

e -(p)+ 7 (k)^e-(p') + 7(k'). 

We will try to follow the general computation that we have developped in the former paragraph. 
However we have to take into account quantum statistics. The general form of the collision term 
(see e.g. ||, H, ||) is 

C[/ 7 ](x,k) = / (-p%)vr + (p) {/ e (p')/ 7 (k') (l + ^) (l " ^) 

-/ e (p)/ 7 (k)(l + ^) [l-^f-))da^\ (82) 

where the factors ^1 + ^y-^ and ^1 — ^ are terms coming from the Bose-Einstein and 
Pauli statistics for the photon and electron respectively p2[| . The differential cross section for 
the Compton scattering is given by 



da kv-*k>p' = I 1 | M fep-^fcV|2 / D kp^k' P '\ 2 (g3) 

4(-p^)' 1 V J ( J 



where \M k P^ k 'p'\ and ^D k P^ k 'p'^j are respectively the matrix element and the two bodies phase 
space element. If we now work in the reference frame of the electron, we have to choose 



— , (84) 



to perform the splitting (|2~l|), from which it follows that the coefficients \M k P~* k ' p '\ 2 and ^D k P~^ k ' p 
are respectively given by 



(D k P- k 'p') 2 = (2vr)V 4 )( P + k- p' -k> + (p> + (k') = \^T 2 ^ (§) (85) 
|M ^V |2 = 16(27r)2m 2 ^ 2 | : + A ^ + + f -2)}g. (86) 
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cjt is the Thomson scattering cross section. The differential cross section is then 

- {l + |A^< + + (87) 

If we use the conjugate process (i.e fcp' — > A/p), we can factorise / e (p)- 
We will also make two following approximations 

• we will take the Thomson limit of the Compton scattering, which implies that 

E E' 
< 1, 

m e m e 

• and we will neglect the quantum statistics. 
With these approximations the collision term can be reset as 

C T [/ 7 ](x,k) =a T J (-p%)7r + (p)/ e (p) (l + \a, v A^ (/ 7 (k') - / 7 (k)) ^, 

from which it follows that the force is 

F"(x») = J C T [/ 7 ](x,k)F7r + (k). (90) 

We can compute this force by using the integrals (28-|30|) given in & [2] and the macroscopic 
quantities ([4^-45), we obtain 

= a T (n efl T^ + p 7 <) . (91) 



Note that this result could have been obtained from the general derivation of the section |5.1| in 
the second approximation (i.e by making tub = m e — > +oo in (|72])). The force on the electron 
is given by 

■^7— »e = — "^e— »7' (^) 

We now need to compute the corrections to this force. They are of two origins, the corrections 
coming from the fact that we do not take the Thomson limit (i.e. corrections in E/m e and 
E'/m e ) and corrections coming from quantum effects (i.e. terms in //). The first one are of the 
same kind that the one for the diffusion (section 5.1). Equation ( |87| ) tells us that ((3 = 0,7 = 
0, a = 3/2) and thus, according to (|77|) , the correction is given by 



= l^ 2 {iV e % - T evX S!? X } (93) 
where we have introduced the macroscopic quantity 

roo r 

I 7 = / / f A E 4 dEdQ, (94) 



m J ft 



We will now deal with the quantum statistics corrections in the Thomson limit. Using the 
general expression of the collision term we have 

F («uant)» = or { /e(p) (/ e ( P ')/ 7 (k) - / 7 (k)/ 7 (k) - / e (pO/ 7 (k'))+/ e (pO/ 7 (kO/ 7 (k)} 

1 + ^A^A') k^ 4 \p + k - p' - k')7r + (p)^ + (p> + (k)^ + (k'). (95) 



To compute this integral, we need to define a whole set of macroscopic quantities related to the 
moments of "/ 2 " . We will not go further here, but we see that such a computation is possible. 
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5.3 Photon-baryon scattering 

If we still work in the Thomson limit, which is equivalent to assume that the electron mass is 
infinite, the collision term coming from the scattering of photons by baryons will be the same 
and thus 

= a T {n w T^ + Pj n;) . (96) 

If we have other fermions in the problem, all interaction between photons and these fermions 
will be described by the same force since we stay in the Thomson limit. The first order term 
will differ from one fermion to the other since it is proportional to m 2 as seen on (p3|). 

5.4 Discussion of the domain of validity 

Three approximations have to be discussed, namely 

• (—p^qn) <C rn 2 AB or A < tuabi i«e the fact that the particles are treated either as non 
relativistic or ultra relativistic, 

• the classical approximation, i.e. the fact that we have neglected the quantum factors 
coming from Bose-Einstein and Fermi-Dirac statistics, 

• the gas approximation 

We have assumed that the particles were non relativistic in the rest frame of the center of 
mass, which ammounts to assuming that the temperature was not too high since we described 
the constituents by massive particles. For a given element this gives a maximum temperature 

n n TUiC 2 

k B 

k B being the Boltzmann constant. The index i emphasizes the fact that the fluids do not have 
to be in thermal equilibrium and can have different temperature (moreover this temperature is 
defined as the statistical temperature). For electrons, we have e * = 6.W 9 K. However, this 
limit can be relaxed if we take into account the corrections in "A/m". For higher temperature, 
electrons can be treated like a radiation fluid, which means that we then used the "ultra- 
relativistic" approximation. 

If we turn to quantum effects, the temperature must not be too low, in order that the medium 
be non degenerate. For fermions of rest mass m and of spin s, this condition (see e.g. P2fl ) is 
satisfied if 

2nh 2 / n \ 2 / 3 
> k B m \2s + 1/ 

In a Friedmann-Lemaitre universe with Oq = 1, the density is equal to the critical density 
p c o ~ 10 -29 h 2 g.cm~ s , where h is related to the Hubble constant via Hq = WOhkm.s" 1 .Mpc^ 1 . 
This leads to an average particles density of no ~ 1 particle. cm~ s for the matter. Thus, for 
electrons (s = 1/2, m = m e ) the former condition and the assumption that n e Q ~ uq gives 

f)> 10 - 16 (-f ■ 

T J \PeoJ 

where To = 2.7K and where we have assumed h ~ 1. Assuming that the electrons behave like 
matter, so that they scale like a 3 , a being the scale factor of the universe, this can be rewritten 
in term of redshift as 

z < 10 16 . 
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For this range of redshifts the gas of photons (or electrons) will behave like a classical perfect 
gas. We see that we have to treat electrons like radiation before we have to take into account 
the quantum effects. The limit for bosons is of the same order of magnitude. 



Since we have computed the general force between uncharged classical particles (section 5.1), 
the formalism can be applied, for instance, to fluid of galaxies in a cluster or of stars in galaxies 
etc.. In such a case one cannot assume that the fluid is composed of particles with the same 
rest mass and each fluid will have a mass spectrum. Since we perform all the integrations on 
T x instead of the mass shell, all the previous results still stand. Note that the approximation of 
non relativistic relative speed is a very good approximation for this kind of fluids. 

The last point we need to discuss concerns the gas approximation, i.e. of the splitting of the 
total fluid (i.e. of energy-momentum tensor) into individual fluids. For that purpose we will 
compare the mean free path, l c say, and the average distance between two particles, d say, for 
the electrons. On the one hand, 

1 1.5 x 10 24 cm 

lc ~ r — \3 — ' 

noT (1 + z) 

where ax is the Thomson scattering cross section (cry ~ 6.65 x 10 -25 cm -2 ). On the other hand, 

— 1/3 

since d ~ n e , using the previous approximate value of n e o, we have 

lcm 

a 



l + z 



The gas approximation is valid if 

l c > d <=^ z < 1.22 x 10 12 . 

Thus, from the time of electron-positron annihilation (T ~ 10 10 -?T,i.e. z ~ 3 x 10 9 ) until the 
time of recombination of hydrogen, it is a very good approximation [||] to treat the content of 
the universe as a nonrelativistic gas plus blackbody electromagnetic radiation. 

Moreover, we know from the theory of the strong interaction that thanks to the "asymptotic 
freedom" , the concept of weakly interacting particles is appropriate for very dense systems [pjj . 
Thus, beyond its conventional range of applicability discussed just above, the hypothesis of 
weakly interacting particles and thus their description by the kinetic theory may be extended 
to early universe. 

As we see on this discussion, it is a crucial point of the theory to assume that kinetic theory 
can be applied and this will have to be checked on any particular case. 



6 Inelastic collisions 

Source terms can arise from many phenomena. For instance if there are some unstable particles 
in the problem, we must take into account decay, fission and fusion. In the early universe matter 
and anti-matter coexisted, which drives us to consider fermion-anti-fermion annihilation. When 
one turns to photons, one major effect has to be considered, namely the Bremstrahlung and the 
recombination. We finish by relating our approach to other work on particles production. 

6.1 General case of fission 

We will try to give the most general term of the source term for fission 

A(p)-B(p') + C(q'). 
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As for the scattering, we start by the kinematics. We work in the rest frame of the decaying 
particle and thus 

u» = —. (97) 
m A 

It is obvious to see that E A = m A an< ^ = 0. This implies that the vector e M = e A is an 
arbitrary unit vector. The vectors p'^ and q'^ are then given by 

p ,fl = \{l + 5 A Bc)p^ + \{l-a A BC + 5ABc) l,2 e^ (98) 
q'^ 1 = \{l-SABc)^-\{l-n A BC + 5ABc) ll2 e^ (99) 

with 5 abc = ( m B ~ m c)/ m A anc ^ a ABC = { m B + m c)/ m A- ^he fission is possible if and only 
if rriA > ins + ma- 
in fact fission is the easiest case of non elastic collision. Since p^ is the only vector of the 
problem and the emission is isotropic in the rest frame of A, it is obvious that 

C[f A ] = C(-p 2 ) = -r A \ (100) 

where ta is a constant representing the lifetime of A. Thus, C[/a] represents the probability of 
decay of A per unit time and the Boltzmann equation with such a collision term (£[/] = —J/t) 
describes the relaxation toward the equilibrium solution f eq = 0. 
The source term is then given by 

e A = I C[f A ]ir + ( P ) = - — ■ (101) 

J TA 

The production of B and C are related to ea by 

e B = e c = -e A - (102) 

The force is 

It can be understood as a "rubbing" coming from the decay. The force on B and C is given by 



TA 



F , = j ^Mp)^ +(p0 , (104 ) 

J TA 

where p' and p are related by the kinematic relations. If we cast ( |98| ) in the equation of F^ B 
and use the fact that the emission is isotropic in the rest frame of A (since is arbitray), we 
get that 

F^ B = ^(l + 5 A Bc)n A /TA and F^ c = l - (1 - 6 abc) n v jJr A . (105) 
It can be checked that, as expected, 

Ft A + F ^B + F Hc = °" 

We must stress that this source term and these forces were computed without any approxi- 
mation. 
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6.2 General case of fusion 

We will try to give the most general term of the source term for fusion 

A(p) + B(q)-C(q'). 

The kinematics of the fusion is analogous to the one of the fission. If we work in the rest 
frame of particle C, we have only to relabel the variables as 

A C, ^ p'», q» -» q^. 
The collision term can be written as 

C[f A ] = ~ J (-^)/A(p)/ B (q)^(p,q,q> + (q)^ + (q'). (106) 



Using the same decomposition that in the section |5.2| , it follows that its general form is 

C[f A ] = -<J J T(mim|,^,e^e^)(-p%)/ 4 (p)/ B (q)^7r + (q). (107) 

Since e'^ is arbitrary, the integration over O' is straightforward and thus 

C[f A ) = -<? J a(mi,m| ) ^)(-^)/A(p)/ B (q)7r + (q). (108) 

We do now the "Thomson" approximation constiting in neglecting (p^q^) compared with Tn 2 A j B . 
In that limit the source term is given by 

e (°) = aa(m 2 A ,m 2 B ) J (~p fl q f ,)f A (p)f B (q)iT + (q)Tr + (p) = (T ABC n B n Afl , (109) 

which means that the reaction rate is proportional to the number of the two reactive constituant 
and to their relative speed. Because of the symmetries we have e A = e B = —tc- Refinements 
can be included if there are any stoechiometric coefficients. 

Let us now turn to the force, 

= aa(m 2 A ,m 2 B ) J (-p%)^/ A (p)/ B (q)7T + (q)7r + (p) = <J ABC n Bv T% . (110) 

If we use the kinematic relations, we can check that 

(1 + 8abc)f!S£ = (1 - 8 ABC )F^, (111) 

which can be understood on two particular cases. If m A = m B then 5 AB c = and F^ A = F^ B , 
i.e. each fluid undergoes the same force because of the symmetry. If m A >ms then S AB c ~ 
and F^ A ~ 0, i.e the variation of impulsion is very small for the fluid A when it merges with 
the fluid B. 

Furthermore, we have 

F^X+F^ = aa(mi,m|)|(-p%)/ A (p)/ B (q)(^ + ^)vr + (q)7r + (p) 

= - f 7Q(mi,m|)|(-p^)/ A (p)/ B (q)p^7r + (p') = -i^ c , (112) 
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because of the kinematics relations and thus 

F (°H^ A + F^g + F®g = 0. (113) 

We can give the correction to the source terms and the force coming from the (—p^q^)- 
dependence in T, which we develop as in the section 5A. This implies a correction to the 
collision terme, C* say, which is given by 

4 1} = / ^WU(p)/B(q)A7T + (p), (114) 

rriAB J 

from which we can compute the correction to the source term, 5eA say, 
(l) _ era 



rriAB 



n BvL K%, whith K% = y/ A (x,p)pVe,vr + (p). (115) 



The correction to the force, F_^ A , is given by 

SF (xg = i p ^ quX f A (p)fB(^ + (p) = -—r A T B ^, (116) 

mAB J rn A B 
where J V A is defined in (f75|). 

6.3 General case of annihilation 

A general annihiliation can be written under the form 

,4(p) + i(q)^ 7 (p')+ 7 (q / ), 

A and A being two fermions. 

An easy way to compute the source term and the force for such a mechanism is to use the 
work we have done on fusion and fission and assume that the annihilation can be seen as 

A(p) + A(q)^B^ 7 (p')+7(q'), 

with the decay time of B being zero (i.e tb — > 0). 
Hence the source terms are 

tA = e A = -e-y/2 

= °AA n A n A^ ( 117 ) 

The force on A and A is the same that the force computed for fusion. The case of the photons 
is a little bit more tedious. A solution will be to compute it from the integral ( p^ ) as we did 
in the former sections. We can however quote that we will have to sum on two photons which 
are travelling in opposite directions and which have the same energy. Thus, by symmetry, the 
total force will be the same as if we had only one particle at the end of the annihilation with 
the 4-momentum 



and with the rest mass 



M 2 = -P 2 . 



This is exactly the situation we have during a fusion and thus, 

F^ A = Ka = vaa^aTa^ F^ = -2a AA n1 A T A)lw . (118) 
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6.4 Photons 

Photons can be produced via Bremstrahlung which can be formally written as 

e(p)^e(p')+7(k). 



If we follow Thorne [15|, the general collision term for such a collision can be written as 

C[f 7 ] = E Xl {E) {G y (E) - / 7 (k)} n e , (119) 

where E = —u^k p and G 7 = %/(£/ 3 % 7 ). r/ 7 and x 7 are the standard emission and absorbtion 
coefficients. If we are in a local thermodynamical equilibrium, then Kirchoff 's law holds and 

G 1 = 2(e E ' T -lY 1 . (120) 



7 

In general, the force on the photons can be written as 

F» = ( — Xl (E)GJE) « + e") n e dEdn 

/rpA 
— (u^ + e^x^f^i^riedEdn. (121) 

This can be computed only if we know the functions Xy{E) and G 1 {E). However we can give 
its general form, which must be 

Ft, = ±- {«(e 7 K - <v(e 7 )} , (122) 

where U(@) and V(0) are two coefficients which depend on the photon temperature, 7 . We 
cannot go further if we want to remain general. However the form has the advantage to be 
covariant and flexible. 

The situation is similar for the source term but, as for the fusion which is in a way very 
similar, we can model it as 

e 7 = ^ and e e = 0, (123) 

T e (fc> 7 ) 

where the rate of emission depends on the temperature. 

If we have other fermions in the problem then the force and the source term for the Brem- 
strahlung induced by these fermions will be alike, with coefficients Uf(@), Vf(Q) and 7y(0 7 ) 



6.5 Recombination 

A special case of interest in Cosmology is the recombination of electrons and protons in hydrogen 
which occurs at the last scattering surface [pTfl . It is however different from the general case of 
fusion since the inverse process is possible, the rate of each process depending on the temperature. 

e(p)+p(q)^^(q')+7(k / )- 

From our study of fusion and fission, it is clear that the source term will take the general form 
(to lowest order), 

£H = e 7 = — e e = — e P 

= a Hl n^n H ^ - a ep n^n p ^ (124) 

The two coefficients, an-y and a ep , are all functions of the temperature. Details on the ionisation 
and recombination processes are needed to have their exact value. Examples of such functions 
can be found in pi, p . 
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6.6 Relation with quantum particle creation 



We want to emphasize that, in our framework, even if there is particle creation via unelastic 
collisions the total energy- momentum tensor is conserved. The situation is then different from 
the one studied in cosmology where the particles production arises from some quantum process 



Such phenomena were phenomenologically described by an effective viscous pressure [ 31 ] and 
a microscopic justification of such an approach was proposed in |32|]. Their work is based on 
the introduction of a source term in the Boltzmann equation beside the usual collision term. 
Unfortunately this source term gives also birth to a force acting on the fluid (they only have 
one fluid). 

If one wants to take into account such quantum creation of particles, this can be achieved 
by using a linear coupling (see Q] section (3-1) 

C[/<] = (-^ + Kx))/*(x,p), 

where u M is, as usual, an arbitrary vector field used to perform the splitting. 

Using the same methods than before it will lead to a source term and a force given by 

f (x) = -^|,(xHx), 
F"(x) = - ^ ' + KxK(x). 



In such a situation the global energy-momentum will not be conserved (see [32] for a detailed 
discussion of this problem). This is beyond the scope of this article. 



6.7 Recapitulation 

Before we give some applications of our formalism, we will sum up the different source terms 
and forces we have computed either in the non-relativistic limit or in the Thomson limit. We 
do not give the corrections than have been computed. In fact, we computed them to show how 
one can give a general form. However, as seen on equations (ff7[), (ff9|), (p3|), ( |115| ) and ( |116| ), 
they imply new macroscopic quantities. 

For any applications, we need as input the cross sections of all the collisions between the 
species we are considering and the lifetime of all the unstable particles (i.e of "cr", "a", 
V, "W" and "V"). 
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7 Application to the theory of cosmological perturbations 



We give here a straightforward application of this general formalism by considering the evolution 
of a gas of photons and electrons interacting via Compton scattering in a Friedmann-Lemaitre 
universe at linear order. The study of the formation and evolution of the cosmic microwave 
background fluctuations begins usually with the study of the radiative transport which is then 
integrated. This has been done by many authors (e.g. 18] f27|| ) in a perturbed Friedmann- 
Lemaitre universe. 

The metric of space-time is given by 



ds 2 = a(rjf (-(1 + 2A)dr] 2 + 2D i Bdx i dr ] + ( 7ij + h^dxUx^ = (g^ + 5g^ v ) dx^dx u , (125) 



if we focus on scalar perturbations. 7^ is the metric of the unperturbed {t = constant}- 
hypersurfaces, D{ is the covariant derivative with respect to 7^, 77 is the conformal time, a the 
scale factor and a prime denotes a derivative with respect to r]. We also introduce h = hij'j 13 . 
Each fluid has an energy-momentum tensor given by 

Ty.v = (P + p)u f ,u u + P(g f , u + 5g^) + (5P + 5p)u fl u u + 5Pg^ + 2(P + p)u (p 5u„ ) +a 2 PU^, (126) 

where P and p are the pressure and density of the fluid in the background and 8p, 5P, 5u^ and 
H^ u are respectively the density perturbation, the pressure perturbation, the velocity pertur- 
bation and the anisotropic stress tensor. Since u^u^ = — 1 then 5uq = and we can write 5v,k 
as 

5u k = a{D k B + v k ). (127) 

We also introduce (? m = dP m /dp m , the sound speed of the matter and 5 = dp/p, the density 
contrast. The law of evolution of the two energy- momentum tensors are, according to equation 

V V T^ = i^U 7 , (128) 
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can be expanded to linear order. Since u^u 7fJi - 
reduces to 

4 / 
F&_^ = -a T n e p y [u: 



1 to first order, the force on the matter ffl 



(130) 



Using the metric ( 125 ) and the above expression of the force, the equation of conservation 
( pljpgD read 



$'m + (! + UJ r, 



Au m + -h' 



1 



1 



(«; + B 1 )' + D % A + ~£>M 7 + jDjTS* 
(v m + fl)' + H(l ~ 3c 2 m )(v m + B)-A 



-aa T n e (vl - 
4 Pl 

-—aa T n e {v^ 

<J Pm 



(131) 
(132) 
(133) 
(134) 



which is the usual result (see e.g. flq| ). Let us stress that since the Compton scattering is 
elastic, the number of photons and electrons are conserved, which can be checked on the first 
two equations. 

If we now turn to the inelastic collision 



A + B ^ C + D 

within the standard model of Cosmology, the equation of conservation for A reads 

V M n^ = -OABn^Bn + cr C Dn^n D ^. (135) 
If we study a small deviation of n A from the equilibrium, then 

V M n^ = -{TABUBeq^ (n A - n^ e J , (136) 

the subscript "eg" meaning at equilibrium. In Friedmann-Lemaitre, this becomes 

n A = -3Hn A - j A (n A - n Aeq ) , (137) 

with j A = dAB^Beq- We recover the general phenomenon of "freezing". If 3H < "f A then the 
system will relax toward the equilibrium, but, if 3H > ja, the system is decoupled since the 
interaction rate is not sufficient to compensate the expansion rate. 

Our formalism, being fully covariant, can be applied to study the same problems in more 
general situations, for instance if we turn to inhomogeneous cosmologies. 



8 Conclusion 

We have developped a fully covariant framework for a system of interacting fluids. We have 
computed general form of the source terms ej, and of the forces Ffii, related to the number flux 
vectors and the energy-momentun tensors T^ v by 

V M < = £ i and VuTr = Ff. 

We have considered the cases of elastic collisions (and the example of the photon-electron scat- 
tering), of inelastic collisions (including fission, fusion, bremstrahlung and recombination) and 
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we have included a possible magnetic force (see table in section 6/7 for a summary of the results 
for all these cases). This computation first required the modelisation of the collision term that 
enters the Boltzmann equation and then the integration of this equation. 

All the quantities have been computed for situations of cosmological relevance (Compton 
scattering, recombination...) and can be used in a wide range of redshift (e.g. < z < 10 12 if 
we just consider photons-electrons and baryons). Since we have the general force for fluids con- 
stituted of massive particles, this formalism is also suited for fluids of stars, galaxies... Moreover, 
thanks to the "asymptotic freedom" , this formalism of interacting gases can be hoped to apply 
to very dense systems, and thus during the early universe. However one has to be careful in 
such an application and must check that the gaz approximation holds. 

This formalism is non perturbative, in the sense that we do not expand the geometry around 
a background spacetime. It is of course perturbative in the sense that we expand and integrate 
the Boltzmann equation to a given order. When expanded to first order (in the perturbation of 
the metric), it reduces to the equation of cosmological linear perturbations. But, our formalism 
can be applied to study the dynamics of a system of gases in a more general context. It will 
be useful, for example, to study the microwave background in inhomogeneous cosmologies or 
compute the sound speed of the photon-baryon fluid during the decoupling. 

Let us stress that we do not assume anything on the fluids but the fact that they are gases, 
i.e that they can be described by a distribution function. When they are perfect fluids, the 
system of equations is closed. Otherwise, we need either an equation of state or an equation of 
evolution (e.g. for II) which can be obtained from higher moments of the Boltzmann equation. 
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